Formation and annihilation of laser light pulse quanta in thermodynamic-like pathway 
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We present a theoretical and experimental study of multiple pulse formation in passively mode- 
locked (PML) lasers. Following a statistical mechanics approach, the study yields a thermodynamic- 
like "phase diagram" with boundaries representing cascaded first order phase transitions. They 
correspond to abrupt creation or annihilation of pulses and a quantized RF power behavior, as system 
parameters (noise and/or pumping levels) are varied, in excellent accordance with the experiments. 
Remarkably, individual pulses carry an almost constant quantum of energy. 
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Lasers operating in a pulse regime can develop more 
than one pulse running in their cavities 1]. That is, at 
any given instant, the optical energy is sharply concen- 
trated around several points along the cavity rather than 
at one point (a single pulse regime) or evenly distributed 
over the cavity (continuous wave operation) . This multi- 
pulse regime can be triggered, for example, by a suitable 
form of absorptive nonlinearity. These lasers exhibit rich 
highly-nonlinear and complicated dynamics, the model- 
ing of which H, 1^ 1^ requires inclusion of nonlinear- 
ities of higher order than the common quartic (Landau- 
Ginzburg-like) term. 

A remarkable feature often exhibited bythese lasers is 
the quantization of the pulse energies 0,13. That is, all 
pulses possess nearly fixed energy, which is almost inde- 
pendent of, for example, the pumping power. These pulse 
quanta can travel with respect to one another, attract or 
repel each other and form ordered structures, such as cou- 
ples or bunches. These phenomena have been l ately re- 
ceiving increasing attention H Ellll El El El El E3 ■ 

Recently it has been shown E3 that the interplay be- 
tween nonlinearities in the laser on one hand and noise, 
such as the inevitable spontaneous emission noise, on 
the other hand, can account for formation and destruc- 
tion of a single pulse in PML lasers. This discontinu- 
ous formation and destruction of pulses was identified 
as a first order phase transition, governed by the bal- 
ance between interaction (due to the nonlinearity) and 
entropy (associated with noise). The study is taken be- 
yond the simple mechanics view of pulses and modes in 
lasers, to a statistical-mechanics picture where a central 
role is played by basic thermodynamic quantities, such 
as "temperature" (noise), entropy and free energy. 

In this Letter we present new theoretical and experi- 
mental results concerning the thermodynamics of pulses 
in lasers. We focus on the very creation and annihi- 
lation of pulses in a multi-pulse regime as the pump- 
ing and the intracavity noise level are varied. We find 
that as noise is increased and/or pumping is decreased, 
pulses successively disappear, one by one, in what is 
shown to be a cascade of first order phase transitions. 
This novel thermodynamic multi-phase light system gives 



rise to a quantized RF power effect. The theory relies 
on an equilibrium-like statistical mechanics model that 
is solved analytically. The qualitative and quantitative 
agreement between theory and experiment is excellent. 

The dynamics of the (classical) electric field in a PML 
laser is commonly described by the complex Ginzburg- 
Landau equation, which can also include high order non- 
linearities 0, especially in the context of multiple pulse 
operation. In our study we wish to include the effect of 
noise as well, which would render the above description 
difficult to handle analytically or numerically. We there- 
fore construct here a much simpler model of the dynam- 
ics, which is capable, as the experiment demonstrates, of 
capturing the key features of noise-dependent formation 
and destruction of pulses. 

We divide the cavity to N equal intervals, such that 
N / T (t is the cavity roundtrip time) is of the order of the 
bandwidth of the laser, which means that the duration 
of a pulse is of order t /N . Then we write an equation of 
motion for the energy Xm at the m-th interval: 



dt 



(1) 



s{x) is the effective nonlinear gain for the energy in the 
interval, t is the long scale time variable over which the 
laser evolves between roundtrips, g is the overall net 
gain (originating from the slow amplifier and effective 
linear losses), Tm is real white Gaussian noise satisfying 
(Fm(t)Fm(t')) = TSmnS{t — t') and the last term is in 
the Stratonovich interpretation [25j . Properly modifying 
g one can always set s(0) = 0. 

When the m-th interval is occupied by a pulse, an 
equation of the form of Eq. Q has been established in 
previous studies and both s{x) and the form of the noise 
term given in Eq. were obtained: The latter for soli- 
tons |20| and the former also for other types of pulses 
|1|. If on the other hand the m-th interval belongs to 
the continuum, its energy is small enough such that non- 
linearities are negligible. Then neglecting s{xm) and the 
dependence of g on Xm in Eq. we are left with the 
same equation of motion that we would obtain for the 
energy of a sample of a band-limited function under the 
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effect of white Gaussian noise. 

We assume that the interaction between different in- 
tervals is weak enough such that the equations of motion 
of Xm-s are not coupled, apart from a global constraint 
of constant total energy V, introduced by choosing g as 
the appropriate Lagrange multiplier 113 • 

While in previous energy rate equation studies |^ the 
energy of the continuum was represented by a single de- 
gree of freedom, in our study this is not appropriate. 
When noise is present the continuum carries essentially 
all the entropy ^^id hence it is crucial that it is rep- 
resented by many degrees of freedom. 

Following the same steps as in Refs. 0,0,^3' 
can show that the invariant measure imposed by Eq. 
is a Gibbs distribution described by the partition function 

Zn{T,V)^ j dxe^ ^^'^''""'^si^x^-V^ (2) 

where S{x) — s{x')dx' and dx denotes integration 
with respect to all x-s from zero to infinity. 

^From Eq. ((2Jl it is clear that the bounded by 

V . If s{x) is an increasing function of x for x < V, which 
can happen when a saturable absorber is present and V is 
small enough, pulses usually do not split. Pulse splitting 
occurs typically when s{x) is an increasing function at 
< X <Xs for some Xg and is a decreasing function for 
X > Xs 13 113, , which is the situation in for example 
in additive pulse mode locking [ij: It is intuitively clear 
that one pulse with a energy much higher than Xs will be 
less favorable than many pulses with a energies of order 
Xs ■ We henceforth assume the above described structure 
of s{x), and show that when V ^ Xg, Eq. predicts for- 
mation of a variable number of pulse "quanta" , i. e. pulses 
with nearly constant energy 22] . Pulse quanta are spon- 
taneously and abruptly created and annihilated when V 
and T are varied. 

The model studied in Ref. is a special case of (2), 
in which the function S is quadratic. There, however, the 
model was derived under somewhat different conditions 
than in the present work. Thus, if s is chosen linear at 
the origin {S is quadratic), the present results recapture 
those of Ref. dSiJ for small V, but we expect that they do 
not describe the actual behavior of the multipulse laser 
at small powers. 

In the thermodynamic theory of our system, the for- 
mation of a single pulse is a first order phase transition 
[itI 0, 0|, and the formation of multiple pulse con- 
figurations is a first order transition between different 
ordered "phases" , reminiscent of structural phase transi- 
tions in solids [2J|. To show the remarkable resemblance 
to the thermodynamic phase picture with the predictive 
power of the statistical-mechanics theory, we refer right 
at the beginning to Fig. ^ The latter shows the theoreti- 
cal "phase diagram" derived below from Eq. Q , and also 



experimental measurements as described below. Each 
phase is labelled by the number of pulses. 
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FIG. 1: Experimental (right) and theoretical (left) phase di- 
agrams. The theoretical graph shows the number of pulses as 
a function of the intracavity energy V and the total noise 
power T, Eq. H12|l . The curves of discontinuity have the 
thermodynamic-like meaning of first order phase transitions. 
The straight lines are asymptotes of the transition lines for 
V>Xs (Eq. CIJ)- 

We proceed to outline the analysis of the statistical 
mechanics model of Eq. which is similar one given in 
detail in Ref. We make the physically appropriate 

assumption that iV 3> 1, i. e. that t is much longer than 
the duration of a single pulse (or equivalently the number 
of modes in the laser, which is of order N , is large). This 
amounts to taking the thermodynamic limit in Eq. 

We start by assuming a specific asymptotic form (as 
N — > 00) for the partition function z''^'^ of configurations 
with m pulses or less 

Z^;p^ (T, P) = (T, 7')e-^(^"(^'^)-i) (3) 

where A^^\t, V) is sub-exponential in N and Fm{T, V), 
the free energy, is the global minimum of 

/™(xi,..,a;„) = -^^S{xj) - lii(v - ^Xj'j (4) 

for nonnegative x values such that "^Xj < V. Let n 
be the smallest m for which the minimum of Fm{T,V) 
with respect to m is attained, and let Xi, ...X„ be the 
corresponding minimizer. Our statement then is that 
Zn approaches Z}^ in the thermodynamic limit, and in 
particular n is the number of the pulses per roundtrip 
and Xi, ...Xn are their energies. If n and the Xj-s are 
to have a finite thermodynamic limit, T has to scale like 
so that T = NT, the total power of noise, has a 
finite limit. The necessity of rendering the parameters of 
the system N dependent is discussed in [l9j. 

Eqs. (|3I4() follow from Eq. ((SJ if one assumes that as 
A'^ — !■ 00 all but a finite number m of the x variables are 
0{1/N). This is justified by a rigorously controlled ap- 
proximation in 1/N based on a recursive equation for Zn, 
as in Ref. . The upshot is that the task of calculating 
thermodynamic properties is reduced to that of finding 
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the minima of / for different values of V and T. In this 
manner, not only the number of pulses per roundtrip is 
obtained, but many other quantities of interest, such as 
the order parameter 



N 



(5) 



which is proportional to the experimentally measurable 
RF power of the photocurrent 26]. 

Since m is finite, finding the minimum of /,« for a spe- 
cific choice of s{x) is straightforward, although in general 
it includes a numerical solution of a set of transcendental 
equations. However for the structure of s{x) we con- 
sider in this Letter, asymptotic expressions for V ^ Xg, 
which is the domain of the multipulse regime, can be 
obtained. Standard minimization techniques show then 
that the minima of are usually obtained at two types 
of configurations: the first has n < m nonzero x values all 
having the same value X , i. e. configurations of n equipo- 
tent pulses, and the second, which is much rarer, is n — 1 
values oix = X and one additional value x = X' < X . A 
sufficient condition for excluding the second type of solu- 
tions is that s{x) increases at least as fast as it decreases: 

a{xi) = s{x2), xi < X2^ \s'ixi)\ > |s'(2;2)|, (6) 

Assuming the first type of solutions, the problem re- 
duces to the minimization of the function 



f{n,x) 



nS{x) 



(7) 



with respect to two variables, n (integer) and x 
(0 < x < V jn). The minimizer X satisfies 



s{X) = 



r 



nX 



(8) 



Clearly X is the common pulse energy and the order 
parameter is Q = nX^ . 

The asymptotic regime relevant to a multi-pulse oper- 
ation, P ^ Xs, is most readily analyzed by considering 
the minimization of the same function / appearing in 
Eq. (O but with n replaced by a real valued variable v. 
Minimizing with respect to x and v together immediately 
gives X = Xt,^ where AT* is the solution of 



S'(A*) = A*s(A*) 



(9) 



A* is the quantized pulse energy. From the assumptions 
made above on six) it follows that A* is unique as long 
as it exists, and A* > Xs (remarkably enough). We use 
the notation s* = s(A*), 5* — S'(A*). Then writing 
A = X^^8 and v = rt+{i/} with n integer and |{t^}| < \, 
and putting back in Eq. JHJ gives 



aI 



s'{X.,){V -vX.,) 



(10) 



if (5 < X*. Since s'(A*) < while s, > 0, indeed 6 < 
A* provided either V — nX ^ A* or n ^ 1. At least 
one of them is certain to hold whenever P ^ A* This 
asymptotic region corresponds to the upper part of Fig.^ 
and is characterized by a nearly quantized pulse energy. 

The phase diagram in the quantized pulse regime is 
very simple: The transition from an [n — l)-pulse config- 
uration to an n-pulse configuration occurs approximately 
when ly is half an odd integer. Using Eq. (jHJ once more 
gives the transition temperature %i{V) 



%{V) = s,V -in-l/2)S, 



(11) 



i. e. the transition curves are approximately equally 
spaced straight lines, as seen in Fig^ 
In order to illustrate our results we chose 



s{x) 



T sin 



2xs 



(12) 



The motivation of this choice is the sinusoidal transmis- 
sivity in additive pulse mode locking [i^l, and although 
s{x), the net gain experienced by apulse, somewhat dif- 
fers from the bare transmissivity [5j, we disregard this 
difference here, since the multiple pulse regime is any- 
way governed by the neighborhood of A*. Results for 
Eq. (|12|) are shown in Fig. ^ for the number of pulses 
and in Figs. Inland 13 for the order parameter Q. 
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FIG. 2: Experimental (right) and theoretical (left, with 
Eq. plots of the order parameter Q (solid line) and the 

mean pulse energy (dashed line) as functions of the noise spec- 
tral power, for a fixed intracavity energy. 

Our experimental study was conducted on a setup sim- 
ilar to a recently reported one |2^ . where the first order 
phase transition associated with the formation and de- 
struction of a single pulse was lucidly demonstrated. It 
consists of a fiber ring laser with PML by nonlinear polar- 
ization rotation technique 0, 0| , with amplified sponta- 
neous emission noise injected from an external source, in 
order to have direct control over the spectral power of the 
additive noise in the cavity. Here we used a shorter laser 
with a roundtrip time of 100 nanoseconds, correspond- 
ing to approximately 20 meters of total cavity length, 
including 4.3-m long erbium-doped fiber amplifier with 
small signal gain of 6dB/m. By proper adjustment of 
the polarization controllers (PCs) PML operation was 
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established with generation of sub-picosecond pulses. As 
observed in a variety of PML fiber lasers P, 0, IE 
excessive pumping (above the self-starting power thresh- 
old) led to the formation of multiple pulses per roundtrip 
that in general were randomly distributed over the cavity. 
However, for certain positions of PCs, stable bunches of 
nearly identical pulses were formed, with approximately 
constant spacing between adjacent pulses (that ranged 
from a few to hundreds of picoseconds, depending on PCs 
positions). 

As the noise or pumping levels were varied, two types 
of responses of the pulse bunch were observed: variations 
in the spacing between adjacent pulses and variations in 
the pulse energy. Therefore, depending on the position 
of the PCs, three distinct regimes of bunched pulse oper- 
ation were obtained. The first and the most common was 
the regime where both types of response were observed. 
In the second regime the multi-pulse bunch contracted 
or expanded while pulse energies remained constant, and 
the third regime was characterized by a fixed bunch pat- 
tern while pulse energies were varied. 

Fig-fflshows the experimental phase diagram measured 
as follows: for several pumping powers the injected noise 
level was raised gradually from zero, the pulses disap- 
peared one by one and the transition "temperatures" and 
average output optical powers were recorded. Such a be- 
havior was previously observed as the pumping power 
was decreased 's^ . The experimental results presented in 
Fig. n were obtained at the first operation regime but the 
structure of the phase diagram was found to be identical 
in all the regimes mentioned above. Fig. ^ demonstrates 
good agreement between theory and experiment. 

Theoretical and experimental plots of the order pa- 
rameter Q (Eq. ©) and the energy per pulse as func- 
tion of the injected noise level (gradually increased from 
zero) for a fixed pumping power are shown in Fig. 13 Ex- 
perimentally they were obtained by measuring the laser 
output with a fast photodiode and an RF power meter 
|26j or a sampling oscilloscope (all having 50GHz band- 
width) correspondingly The pulse energy is nearly con- 
stant, with deviations of about 5%. These deviations are 
well described by the theory. The results of Fig. [3 were 
obtained at the third operation regime (constant spacing 
regime described above, where pulse energy changes are 
most pronounced). 

Fig. Ul shows additional theoretical and experimental 
plots of the order parameter dependence on noise for both 
increasing and decreasing of the noise level. Typically to 
first order phase transitions, the system exhibited hys- 
teresis: The number of pulses at any point {T,V) de- 
pends on the precise path that led to it. In particular, 
increasing T leads to a different Q{T) curve than decreas- 
ing it, as seen in Fig. |3 The hysteresis can be theoreti- 
cally estimated using an analog to the Arrhenius formula: 
the system dwells in a meta-stable phase, corresponding 
to a local minimum of /, as long as this minimum is 
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FIG. 3: Experimental (right) plot of the order parameter Q 
as a function of increasing (solid line) and decreasing (bro- 
ken lines) noise spectral power, demonstrating three different 
hysteresis paths, with excellent agreement to the theoretical 
plot (left), with the choice s(x) = aix — a^x'^ + a^x^. A mo- 
tion picture demonstration of the experimental results can be 
views in EPAPS document no [XXXX]. 

"deep" . Since the lifetime of a metastable phase is ex- 
ponential in the barrier surrounding the corresponding 
local minimum of /, it can be much longer than the time 
over which the system parameters are varied. 

We conclude by noting how powerful the combination 
of statistical-mechanics and laser physics can be, leading 
us to a new view and findings that can be significant to 
both fields. 
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